Charge Density Bounds in Superconducting States 
of Strongly Correlated Systems 
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Charge density bounds of knotted and linked vortex states in two-component Ginzburg-Landau 
model are considered. When the mutual linking number of vector order parameter vortex lines is 
less than the Hopf invariant, these states have the lower-lying energies. It is shown that a set of 
local minima of free energy contains new classes of universality which exist in a hole density range 
limited at both finite ends. 
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Introduction. - A tangle of vortex filaments is a system 
which attracts attention due to several reasons. Along 
with the coherent state, which is the background of this 
vortex field distributions, the filament system also con- 
tains a disorder combination due to free motion of its 
fragments and a topological order because of the effects 
of knotting and linking of its separate parts. The degree 
of this order can be characterized by the linking number 
of vortices. Persistent linked field configurations along 
with a large value of interaction energy present the basic 
content of what we mean when talking about the origin 
of strong correlations. 

The study of soft condensed matter, whose universal 
behavior is determined by topological characteristics, is 
an active area of research in physics [Q, [|, ||, Q| and be- 
yond H . Indeed, the dynamics of entangled vortex states 
is recognized as one of the most challenging problems 
of modern condensed matter physics. The emergence of 
novel analytical methods ||, || to treat nonlinear field 
equations create now an unique opportunity to under- 
stand the dynamics of entangled vortex states. 

The vortex dynamics can be understood in detail in 
the framework of basic models, such as the Ginzburg- 
Landau functional or more complicated multi-component 
systems. Main theoretical approaches for analyzing uni- 
versal behavior in these models are based on the methods 
of the topological field theory || and nonlinear dynam- 
ics j?J . The universality classes are defined both by sym- 
metry and topological characteristics of the background 
nonlinear fields, while a system with infinite number of 
degrees of freedom contains in the ground state only the 
finite number of symmetrically invariant physical states. 

The aim of this paper is to elucidate the physical mech- 
anism leading to formation of the charge density range 
bounds for superconducting states with a set of num- 
bers determining the knotting and linking degree of the 
fields that take part in the description of the coherent 
state. The use of the Ginzburg-Landau functional and 
the methods of the topological gauge field theory gives 
us the hope that the obtained answers are universal. A 
combination of novel methods will make it possible to 
advance in description of the processes which are deter- 
mined by the contribution of topological excitations. 



Model. - We will use the Ginzburg-Landau free energy 



d 3 x 



E- 

^ 2m 



2e 

hd k + i—A k * Q 



E(-m* q i 2 + > q i 4 ) + C 



with a two-component order parameter 
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satisfying the CP 1 condition, |xi| 2 + |x2 1 2 = L This 
model is used in the context of the two-gap supercon- 
ductivity || and in the non-Abclian field theory (|, [To| . 

It has been shown in paper j^] that there exists an 
exact mapping of the model (1), (2) into the following 
version of n-field model: 



F= d 



^(d k nf + (d k p) 2 + ±p^ 



+ (F ik - H ik f + V{p,n 3 ) 
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To write down Eq. (3), dimensionless units and gauge 
invariant order parameter fields of the unit vector n = 
X&x, and the velocity c = J/p 2 have been used. Here 
X = (.X11X2) an d c - Pauli matrices. The effective cou- 
pling constant 512e 2 /fi.c in this case Q has the order of 
unity. The full current J = 2p 2 (j — 4A) has a paramag- 
netic (j = i[xiVxi — c.c. + (1 — > 2)]) and a diamagnetic 
(— 4A) parts. Besides, in Eq. (3) Fi k = d L c k — d k Ci, 
H ik = n- [diii x d k n] := dia k - 9 fe a 4 . 

Setting in Eq. (3) c = we come back to the model of 
pd[ . The numerical study of the knotted configurations 
of n-field in this model has been done in |l^, [l3| Q . The 
lower energy bound in this case 

F^32tt 2 |Q| 3 / 4 (4) 

is determined jl5|, [l6|, [TtJ by the Hopf invariant, 
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At compactification R 3 — > S 3 and neS 2 , the integer 
Q S 7r 3 (5 2 ) = 2 shows the degree of linking or knotting 
of filamental manifolds M E S 3 , where the vector field 
n(ir, y, z) is defined. In particular, for two linked rings 
(Hopf linking) Q=l, for the trefoil knot Q=6 and etc. 
Significant point Jl8| is as follows: tt 3 (CP m ) = at M > 
1 and ^(CP 1 ) = ~ks(S 2 ) — 7L. In the latter case the 
order parameter (2) is two-component one ^ and linked 
or knotted soliton configurations are labeled by the Hopf 
invariant (5). 

Energy bounds. - Let us assume that p = po can be 
find from the minimal value V(po) of the potential V(p), 
but the velocity c does not equal zero. Equation (3) in 
this case has the following form: 



F — F n + F c — F int 
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It is seen from Eq. (6) that a superconducting state 
with | cj <C 1 has the energy which is less than the mini- 
mum in Eq. (4). To find the lower free energy bound in 
the superconducting state with c =/= 0, it has been proved 
|fH| that the following inequality for free energy parts 
takes place 



Fl^Fl' 2 > (32^ 2 ) 4 / 3 |L| 



where 
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is the degree of mutual linking |2lJ] of the velocity c lines 
and of the magnetic field H = [V x a] lines. Quantities 
like (5), (8) arise as first integrals in the theory of an 
ideal or barotropic fluid pl[ . In general, the integrals in 
Eqs. (5), (8) could be defined by the asymptotic linking 
numbers p2|| . 

It follows from the Schwartz-Cauchy-Bunyakovsky in- 

1/2 1 /2 

equality that 2P„ • P c ^ F int . Setting the boundary 



value Fmt into Eq.(6), we get F„ 



(F 7 ] /2 - P c 1/2 ) 2 . 

The mimimum value 32tt 2 Q 3 / 4 of the function F n and 
Eq. (7) lead El to 



F^32^ 2 |Q| 3 / 4 (1-|L|/|Q|) 2 



(9) 



The trivial case Q = should be considered after the 
limit L — 0. Let us pay also an attention to the self-dual 
relation F n = F c which follows from F min . 
In general, the Hopf-Chcrn-Simons matrix, 
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can be defined in our case by linking numbers Q 6 Z 
of spin and spin-charge ({L, Q'} ^ Z at all) degrees of 



freedom. In Eq. (10), K a p is a symmetrical matrix (JJ 



L) 



and a} 



c,-. Following the approach of 



|19[ which is based on employ of the Holder inequality 
chain as well as on the Ladyzhenskaya pof inequality, 
one can find that 



P„ 1/2 P C 5/6 ^ (167r 2 ) 4 / 3 |L'| 



(11) 



The distiction of the coefficient in Eq. (11) from Eq. (7) 
results from the charge 2e of pairs (due to the coefficient 
1/4 in the firt term of the free energy part F c ). 

It follows from Eq. (9) that for all numbers L < Q 
the energy of the ground state is less than that in the 
model described in [JlTJ , for which the inequality (4) is 
valid. The origin of the energy decrease can be easily 
understood. Even under the conditions of the existence 
of the paramagnetic part j of the current J, the diamag- 
netic interaction in the superconducting state consumes 
its own energy of the current and a part of the energy 
relating to the n-field dynamics for all state classes with 

Charge density. In the context of the high- 

temperature superconductivity problems |2^, |24|, EM |6| , 
Eq. (2) has the sense of the condition of the correlated or 
linked factorization of charge and spin degrees of freedom 
(from the gauge invariant order parameter point of view 
(3)). In fact, the constant value of the charge density 
plays the role of a tuning parameter of the system. In 
this case, the order parameter n describes the distribu- 
tion of spin degrees of freedom, while the order parameter 
c contains the contribution both of spin degrees of free- 
dom Xa in the paramagnetic part of the current and of 
the U(l) charge degrees of freedom in the diamagnetic 
part of the current J. Due to such mixing, the gauge 
invariant order parameters n and c differ principally and 
topologically. If the mutual linking index is defined by 
the current J, the measure in Eq. (8) depends on the 
charge density. The vector c = 3 / p 2 is normalized to the 
charge density and, as distinct from n, belongs to the 
noncompact manifold. This leads to the fact that the 
Hopf numbers for it are not integer valued ((L, Q') ^ 2), 
and the linking numbers arc included into the factor of 
Eq. (9) in the form of the ratio L/Q. 

Let the parameter po change in some range. Since all 
terms in Eq. (3) are of the same order, the velocity c and 
therefore the number L decrease as po increases. In this 
case, if L is rather small, the smallest superconducting 
gap goes down with the growth of Q at the background 
of large value 32tt 2 Q 3 / 4 of the spin gap. 

As po decreases, the following effect takes place. The 
radius R of the compactification R 3 — > S 3 being propor- 
tional to Pq , increases till it exceeds a certain critical 
value, R cr . At R > R cr the Hopf map |2?J is unsta- 
ble relative to infinitesimal perturbation of vortex linked 
field distribuions. As a result, the U(2) symmetry asso- 
ciated with identical Hopf map is spontaneously broken. 
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This means that the topological solitons, instead of be- 
ing spread out over the whole of S 3 , are then localized 
around a particular point (the base point of the stereo- 
graphic projection) and collapse into localized structures 
. This picture corresponds to the existence of an op- 
timal value of the relation L/Q. Under such a key condi- 
tion of the restoration of compactness of the base mani- 
fold, inhomogeneous superconducting state appears. As 
the parameter p decreases further, the free energy F 
due to the inhomogeneity term (dkp) 2 increases again 
and suppresses the superconducting gap. 

Up to now the vector A characterized the internal 
charge U(X) gauge symmetry. If we apply an external 
electromagnetic field the vector potential A equals the 
sum of the internal and the external gauge potentials. 
As a result, due to the diamagnetism of the supercon- 
ducting state, the velocity c decreases. Like in the case 
of increasing p , this leads to suppression of the super- 
conducting gap. The answer to the question on the exis- 
tence of full or partial Meisner screening in these states 
depends on the result of the competition of contributions 
from neutral j and charged — 4 A parts to the full current 
J. 

The soft case p ^ const both for c = [^9| and c ^ 
arouses certain interest. It is more complicated due to 
some reasons and will be considered in a separate paper. 
The point y/F c /F n — 1/4 where the fields of the in- 
equalities (7), (11) intersect, deserve an additional care- 
ful study as well. The equality in Eq. (9) under this 
remarks should be understood as an ideal limit depend- 
ing on topologial characteristics of knots and links only. 

Discussion. - In the (3+0)D case of the free energy (3), 
Hopf invariant (5) is analogous to the Chern-Simons ac- 
tion (k/An) J dtdrxEnvxa^dvOLx determining strong cor- 
relations of (2 + l)D modes @ |l| at k ~ 2. In pla- 
nar systems, this coefficient has the sense of braiding of 
the excitation world lines. In particular, for the semion 
k = 2. Keeping in mind the relation of spatial dimen- 
sionality of the systems in their quantum and statistical 
descriptions, we note that the (2 + 1)D dynamical case 
k = 2 of the open world line ends of excitations is equiv- 
alent to the compact (3 + 0)D statistical example of the 
Hopf linking Q = 1 in Eq. (3) (see Fig. |l|). 

Similarly to the description of the states in the frac- 
tional quantum Hall effect with the filling factor v = p/q 
and p,q £ Z, the Eq. (9) depends on L/Q. The Hopf 
invariant Q € Z enumerating gauge vakua [^8| , is equiv- 
alent to the degree q of the ground state degeneracy and 
L plays the role of the filling degree p of the incompress- 
ible charged fluid state. We have stressed already that 
the number L is, generally, not integer valued. From 
this point of view the multiplier (1 — L/Q) in Eq. (9) is 
equivalent to the Hall filling factor 1 — v for holes. (Let us 
remark that the mentioned parameter has the sense 
of the particle density). The distinction of our system 
from the fractional quantum Hall effect state is that due 




FIG. 1: Examples of n- and c-field line configurations in 
the form of Borromean rings with vanishing pairwise linking 
numbers K a p (on the l.h.s.) and in the form of the closed 
chain of correlated charge-spin semion loops with Q = 1 (on 
the r.h.s.). 

to compressibility of the superconducting state, where 
the charge gauge f (1) symmetry is broken, the effective 
number L of the charge degrees of freedom is continuous. 
It plays the role of an interpolation parameter ]32] | which 
connects sectors with L ^ Z and I € Z in one and the 
same universality class of Pfaffian states. Apparently it 
is true for all Pfaffian- like states with Q > 1 and numbers 
L < Q. 

In conclusion, we have studied the charge density 
bounds of the superconducting states using the CP 1 ver- 
sion of Ginzburg-Landau model under the conditions of 
the existence of linking and knotting phenomena of the n- 
and c-fields being the gauge invariant order parameters 
of the considered system. We have shown that the super- 
conducting states exist in the finite range p\ < p§ < p\ of 
the charge density values. In this range, the field config- 
urations with field lines characterized by semion values 
of the linking numbers are preferable. 

I would like to thank A.G. Abanov, E.A. Kuznetsov 
and G.E. Volovik for advices, V.F. Gantmakher for the 
crucial remark, and L.D. Faddeev, G.M. Fraiman, A.G. 
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was supported in part by the RFBR under the grant No. 
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